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ABSTRACT 


Shock trials for naval vessels are a requirement for each new class of surface 
ships in the U.S. Navy. With understanding the technology of underwater shock analysis 
and considering the rising costs of conducting actual shock tests, computer simulation of 
shock trials is becoming more and more attractive. Unfortunately, finite element models 
can be quite large and require sufficient amounts of computer memory and time to run a 
shock analysis. This thesis investigates the effects of reducing the element size of a ship- 
like box model subject to an underwater explosion. Known as smearing, this process 
combines the density and stiffness properties of the removed elements into the remaining 
material of the model. Positive results from computer simulation could greatly affect the 


manner in which shock trials are conducted with future ship classes. 
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I. INTRODUCTION 


A BACKGROUND 

During the Second World War, the United States Navy began to investigate more 
thoroughly the damage sustained by surface ships due to impacting shock waves, or 
pressure pulses, produced from underwater explosions such as those of non-contact 
mines. These pressure pulses, applied to a large area of the ship’s hull, can have a serious 
negative effect on the ship structure as well as the equipment located onboard. The 
purpose of these tests was to improve the survivability of, while minimizing damage to, 
the ship’s hull and onboard equipment and systems. Since then, the Navy has developed 
guidelines and specifications for shock testing new classes of ships and the equipment 
located in them. Known as shock trials, these tests are required for all surface ship 
classes that must be capable of operating in a combat environment. Shock hardening 
criteria apply to onboard systems and equipment that perform or directly support mission- 
essential functions on these ships. The goal of shock hardness testing is the development 
and installation of system improvements to prevent or minimize the broad degrading 
effects of shock on mission performance [Ref. 1]. Further guidelines on shock hardening 


of surface ships are outlined in OPNAVINST 9072.2 (reference 1). 


Unfortunately, conducting ship shock trials can be time consuming and expensive. 
Initial planning stages for the shock trial of the USS John Paul Jones (DDG 53) began 
four years prior to the actual test date while the ship was still undergoing construction at 
Bath Iron Works, Maine. The entire endeavor involved over 50 government agencies. 
Originally scheduled for February 1994 when weather conditions are most favorable, the 
shock trial was delayed 3 months due to a lawsuit filed against the Navy by 
environmentalist groups concerned over the well being of sea life in the testing area. 
When testing occurred in June 1994, only two of the four required tests could be carried 
out because of inclement weather. The remaining two shock tests could not be performed. 
Further modifications to the ship’s schedule to accommodate the two tests were not 
feasible since the three-month delay had already affected the ship’s post trial delivery 


date and deployment preparations [Ref. 2]. 


With the advent and ongoing advances in computer technology, finite element 
modeling has become more and more refined. In view of the time and cost involved in 
planning and conducting an actual shock trial, simulating these tests on a finite element 
computer model would be desirable in many aspects. Finite element modeling codes 
such as TrueGrid®, MSC/PATRAN, and I-DEAS have enabled generation of detailed 
finite element models in a timely manner. In analyzing finite element models, higher 
detailed model meshes provide more reliable results than do less refined meshes. Highly 
detailed models, of course, do require valuable computer memory and time when 
performing structural analysis. Reducing the complexity of a finite element model mesh 
is one method in which to solve the problem of time consuming and expensive computer 
usage. Even though a less refined mesh may affect the results of the analysis as 
mentioned above, acceptable results can still be obtained depending on how much the 


model complexity is reduced. 


B. SCOPE OF RESEARCH 

This paper investigates the effects of reducing the complexity, or element size, of 
a finite element mode] subject to an underwater explosion. The reducing process is 
known as smearing and combines the density and stiffness properties of the removed 
elements into the remaining material of the model. The model to be considered in this 
study is a ship-like box model with two bulkheads, a keel, and beam stiffeners. Analysis 
of the model response is conducted using the Underwater Shock Analysis (USA)/LS- 
DYNA. The purpose of this thesis is to compare responses of smeared models to those of 
the original in order to ascertain the reliability of this model reduction method. With 
reliable results, computer simulation of ship shock trials could be a dependable, cost 
effective, and time efficient manner for validating surface ship shock hardening 


requirements. 


Il. UNDERWATER EXPLOSIONS 


A. UNDERWATER SHOCK PHENOMENA 

The sequence of events for an underwater explosion is fairly complex. When a 
high explosive (TNT, RDX, HBX-1, etc.) 1s detonated, the original solid material of the 
explosive is converted into a high temperature and pressure gas within nanoseconds. The 
pressure wave that is formed originates in one section of the explosive and propagates 
throughout the remainder of the explosive. As this pressure wave propagates, it initiates 
the chemical reaction that creates more pressure waves. The pressure wave velocity 
steadily increases within the solid explosive until 1t exceeds the speed of sound in the 
explosive. Therefore, a shock wave is produced. The shock wave propagates through the 
solid at a constant speed and then, with the high temperature and pressure behind the 


shock front, into the surrounding medium [Ref. 3]. 


The high-pressure gas that results from the explosion rapidly expands outward in 


a radial manner, Fig. 1, and imparts an outward velocity on the surrounding water as 






High Pressure Gas Bubble 


Charge Diameter 


Shock Wave Front 


Figure 1. Gas Bubble and Shock Wave from an Underwater Explosion 


well. Initially, the pressure 1s much greater than the atmospheric and hydrostatic pressure 
that opposes it and is therefore compressive in nature. At detonation, the pressure rise is 
discontinuous and decays exponentially with time as shown in Fig. 2. Duration of the 


pressure disturbance lasts only for a few milliseconds. The shock wave is characterized 


34,000 psi 


3,400 psi 
t 2,200 psi 
ue 4 se PS! 460 psi 
5! 45' 50’ 495' 500’ 


Figure 2. Shock Wave Profiles From a 300 lb. TNT Charge [Ref. 3] 


by a propagation velocity that is several times that of acoustic velocity, 5,000 ft/sec, near 
the charge which then falls rapidly to acoustic velocity as it travels outward in the water. 
Additionally, the pressure profile of the shock wave is proportional to the inverse of the 
distance from the charge, 1/d, and the wave profile gradually broadens as it spreads out 
[Ref. 3]. Empirical equations have been determined to define the profile of the shock 


wave: 


t=) 


as ae 8 Gs) 





pe aes) 


0=KW: _" (msec) 


where P is the magnitude of the pressure wave at time t, Pmax is the peak magnitude of the 
pressure in the shock front, t-t; is the elapsed time after the arrival of the shock, @ is the 
decay constant which describes the exponential decay of the pressure wave at one-third 
of its maximum value, R is the distance from the explosive to the target in feet, and Aisa 


constant that depends on the charge type. 


Depending on the charge location relative to the surface and the bottom, other 


effects are characteristic of an underwater shock. Surface cut-off is a tension wave, as 


opposed to a compressive wave, that is produced from the rarefaction of the shock wave 
off of the surface. Bottom reflection is the bouncing of the shock wave off of the bottom 
of the body of water; a compressive wave. Refraction causes the shock wave to travel 
through the bottom of the body of water before emerging again; also a compressive wave. 
Surface effects also occur as a result of an underwater explosion. A spray dome appears 
first which is caused by the shock wave. Water plumes follow the spray dome and are a 
result of the bubble pressure pulses. A bubble pressure pulse is produced when the gas 


bubble collapses during its oscillation and migration to the surface [Ref. 4]. 


B. FLUID-STRUCTURE INTERACTION 

When an object such as a ship or submarine is in the vicinity of an underwater 
explosion, the shock pressure pulses produced by the explosion impinge upon the surface 
of the structure. A fluid-structure interaction takes place as the pressure pulse acts upon 
the flexible surface of the structure. The entire structure is excited by the pressure pulse 
and responds in a dynamic manner. Taylor flat plate theory is the simplest case of fluid- 
structure interaction, first introduced by G. I. Taylor. In this case, an infinite, air-backed 


plate is acted upon by an incident plane shock wave as shown below in Fig. 3. 


P ie 
Air \ 1 Air 2 
BackedN Backed 
4 2 
Shock Wave Approaching Shock Wave Reflecting 


Figure 3. Taylor Flat Plate Theory: Incident and Reflected Pressure Waves 


Once the shock wave strikes the plate, a reflected shock wave leaves the plate. 
According to Newton’s second law of motion as shown in Eq. (2.1), 


5 


Tees ea 2.1) 


dt 


where m is the mass of the plate per unit area, u is the velocity of the plate after being 
subjected to the shock wave, P(t) is the incident wave pressure and P>(t) is the reflected, 
or scattered, wave pressure. Define the fluid particle velocities behind the incident and 
reflected shock waves as uj(t) and u2(t). The velocity of the plate is then defined by Eq. 
(2.2), 


u(t) =u, (t)—u,(t) (2.2) 


For the one dimensional plane wave, the wave equation is P = pCu. It follows that the 


incident and reflected shock wave pressures, Eqs. (2.3) and (2.4), are 
P =pCu, (23) 
P, =pCu, (2.4) 


where p is the fluid density and C is the acoustic velocity in water. Substituting the 
above pressure Eqs (2.3) and (2.4) into the velocity Eq. (2.2) results in the incident shock 


pressure being defined as Eq. (2.5) 
P(t)=P.,,e ° (2.5) 
Solving for the reflected shock pressure yields Eq. (2.6), 
P,(t)=P —pCu=P__e ®*—pCu (2.6) 


where t is the time after the shock wave arrives at the target. Now the equation of 


motion, Eq. (2.1) can be rewritten as Eq. (2.7), 
du = a 
on +pCu=2P _e (25) 
t 


which is a first order, linear differential equation. The solution, u(t), of the differential 


equation is expressed in Eq. (2.8) as 


2Bt = 
oP a le —¢ F (2.8) 


y=-e___ 

m(1—B) 
with B = pC@/m and t>0. The total pressure that impinges on the plate is defined as Eq. 
(2.9), 


——e °] (2.9) 


As the value of B becomes larger, as in the case of a lightweight plate, the total pressure 
will become negative at a very early time. However, since water cannot support tension, 
negative pressure cannot exist. Therefore, as the water pressure reduces to vapor 
pressure at the surface of the plate, cavitation occurs. At this point, the pressure in front 


of the plate has been cut off and the plate has reached its maximum velocity [Ref. 4]. 


The total pressure relationship shown above in Eq. (2.9) is a critical component in 
solving the fluid-structure interaction problem. In 1971, Thomas L. Gears described a 
matrix of differential equations in time for the approximate treatment of acoustic fluid- 
structure interaction. It is called the Doubly Asymptotic Approximation (DAA) and is 
accurate at both low and high frequencies and at early and late times. The DAA 
represents the surrounding fluid of the structure through the interaction of state variables 
pertaining only to the structure’s wet surface [Ref. 5]. The fluid equation, Eq. (2.10), of 
the DAA is 


[M,]{P.}+ pclA, ]{p,} =eclM, iu, } (2.10) 


where [Mg] is the symmetric fluid mass matrix for the wet-surface fluid mesh, {p2} and 
{p,} are the nodal pressure vector and its first time derivative of the scattered wave, c is 
the acoustic velocity of water, [A¢] is the diagonal area matrix associated with the fluid 
elements, and {u,} is the scattered wave velocity vector. An object in the fluid will have 


a structural response defined as follows in Eq. (2.11), 


[M.]{X}+[C,] {x} + D8, ]{x} = tf} (2.11) 


where [Msg] is the mass matrix, [C;] 1s the dampening matrix, [K,] is the matrix, stiffness 


{x} is the acceleration vector, {x} is the velocity vector, and {x} is the displacement 


vector of the structure and {f} is the external force vector. Compatibility at the wet 
surface is described by Eqs. (2.12) and (2.13) 


(f} = -[G][A, ]({P,$ + (Pt) (2.12) 


and 


[G]’ {x} = {u,} + {u,} (2.13) 


where [G] is the transformation matrix that relates the surface nodal forces of the fluid 
and structure. The “T” superscript indicates the transpose of the matrix. Substituting 
Eq. (2.12) into (2.11) and solving for {uz} in Eq. (2.13), taking the first time derivative, 
and substituting into Eq. (2.10) results in Eqs. (2.14) and (2.15), 


[M,]{X} +[©, ]{x} +[K, ]{x} =-[G]IA, JCP, } + {P.}) (2.14) 
and 


[M,]{p.}+pc[A,]{p,} =peM, ([G]’ {x} — {4,}) (2.15) 


The USA code solves Eqs. (2.14) and (2.15) simultaneously by using a staggered solution 
procedure that is unconditionally stable with respect to the time step used [Ref. 6]. Once 
this system of equations is solved, desired response results such as displacement, 


velocity, and acceleration can be studied. 


Il. MODEL CONSTRUCTION, TESTING AND SMEARING 


A. MODEL CONSTRUCTION AND TESTING PROCEDURE 
Model construction and testing involves pre-processing, conversion, analysis, and post- 
processing programs. A flow chart of the model building and testing procedure is shown 


below in Fig. 4. 
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Figure 4. Flow Chart. Model Construction and Testing 


The box-like ship model for this study was created using a finite element mesh 
generator known as TrueGrid®. This computer program produces hexahedral and 
quadrilateral (quad) elements from geometric model data that is input by the user [Ref. 
7]. After necessary material properties are defined, TrueGrid® writes the model to an 
output file in NASTRAN format. The NASTRAN input file, or deck, is then entered into 
PATRAN, a finite element mesh generator, analyzer, and pre- and post-processor [Ref. 
8]. In PATRAN, the NASTRAN deck is then set up for normal modal analysis in order 
to verify a correct dynamic response of the box structure to ensure that the model is 
constructed properly. Modal analysis is performed by NASTRAN, a finite element 
analysis code that allows for a number of different analysis types (Ref. 8]. The next step 
is to convert the NASTRAN deck into a LS-DYNA input deck. LS-DYNA is a finite 
element code that provides the means for analyzing structures under a variety of load 
conditions [Ref. 9]. For this particular study, LS-DYNA3D was coupled with the USA 
code to obtain the results to the underwater shock phenomenon. The USA code solves 
for the transient response of a partially or totally submerged structure subjected to an 
acoustic shock wave [Ref. 6]. Once the transient analysis is complete, the response of the 
structure is converted into ASCII files using LS-TAURUS. This program is used for 
post-processing two and three dimensional finite element analysis codes [Ref. 9]. 
Responses from the structure were then graphed using the Underwater Explosion 
Research Division (UERD) Tool, a PC based program for displaying a variety of shock 


responses. 


The box structure analyzed in this research, shown in Fig. 5, contained two 
bulkheads and a keel. It was 120-in long, 24-in wide, and 24-in deep. The bottom, sides, 
and bulkheads of the box were made of 1/4-in steel plating with a weight density of 0.284 
lbf/in? and an elastic modulus (Young’s modulus, E) of 30 x 10° psi. Steel stiffeners with 
the same mass density and elastic modulus as the plate material were added to the 
bottom, sides, and bulkheads to increase the model rigidity. Their cross section 
dimensions were 0.125-in thick by 2-in high. Keel cross section dimensions were 0.25-in 
thick by 6-in high. Four lumped masses (0.138-lbf s*/in* each) were added to the keel in 


order to make the draft of the box structure one-foot. The model’s finite element mesh 
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Figure 5. Box Model Dimensions. A) Bottom B) Sides C) Ends/Bulkheads 
D) Stiffener and Keel Cross Sections E) Model Cross Sections 
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contained quad elements and beam elements of varying size depending on their location 
in the mode]. Quad elements on the two ends and bulkheads were 4-in by 4-in and beam 
elements were 4-in long. The bottom and sides of the two end compartments contained 
4-in by 9-in quad elements with 4-in beam elements in the lateral direction and 9-in beam 
elements in the longitudinal direction. The sides and bottom of the center compartment 
consisted of 4-in by 9.6-in quad elements with 4-in beam elements in the lateral direction 
and 9.6-in elements in the longitudinal direction. In all, the box structure contained 386 
nodes, 378 quad elements, 615 beam elements, and 4 point elements that were necessary 
for the added lumped masses. The structure finite element mesh is shown in Fig. 6 and 


stiffener placement is displayed in Fig. 7. 


B. SMEARING 

As mentioned in the first chapter, the complexity of the box model was reduced 
by a process known as smearing. The density and stiffness properties of the removed 
elements were combined into the remaining material. In the case of a stiffened flat plate, 
smearing combines the elastic modulus and density of each stiffener with those of the 
base plate to form an equivalent orthotropic base plate. Physically, the equivalent, or 


smeared, plate will have the dimensions of the original base plate, Fig. 8. 


Given a stiffened, flat plate with dimensions as shown in Fig. 9, the flexural 


rigidity of the equivalent plate is defined in Eqs. (3.1) and (3.2) as [Ref. 10], 


(El), = (ED, + ¥. (ED. 3.1) 


i=] 


OT 
EI, =E,I, + ¥ (EI, /a) (3.2) 
i=) 


where (EI), (ED, and (ED), are the flexural rigidities of the equivalent plate, the original 


base plate, and the stiffeners, respectively, and n is the number of beams smeared into the 
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Figure 7. Stiffener Placement for the Box Structure 








plate in the direction of interest. The moment of inertia of the equivalent and original 
basesplates 1issia—sipaln — ia Geae and the moment of inertia of a stiffener is I, = 
(oe fie) + b,h,d*. It then follows that the equivalent modulus of elasticity for the 


smeared plate is as shown in Eq. (3.3) 


© 12E 1 (1-v’) 


E,=E, + : 
i=l a(h, ) 


€ 


(3:3) 


where E, is the modulus of elasticity of the equivalent plate, E, is the modulus of 
elasticity of the orginal base plate, E, is the modulus of elasticity for an individual 
stiffener, I, is the moment of inertia for an individual stiffener about the centroidal axis of 
the base plate, v is Poisson’s ratio, a is the distance between stiffeners, and hy is the 
thickness of the base plate. However, in this study, it was assumed that the equivalent 
plate thickness, he, is the same as the original plate thickness, hy, and I, = I, = wh,°/ 12 
where w is the width of the plate being smeared. The equivalent modulus of elasticity for 


the smeared plate used in this study is then defined by Eq. (3.4) 


E, =E, +e (3.4) 
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Figure 8. Smearing of a Stiffened, Flat Plate 
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The density for the equivalent plate is defined by Eqs. (3.5) and (3.6) as follows [Ref.10], 





DAY pee Pidin Yp.v (3.5) 
and 
OV 
P. =P, PY : (3.6) 


where P¢ is the density of the equivalent plate, pp is the density of the base plate, p, is the 
density of an individual stiffener, vp is the volume of the base plate, and v, is the volume 


of an individual stiffener in the direction being smeared. 





Figure 9. Smeared Plate Variables 


The box model was reduced twice in order to study the response of the smearing 
method mentioned above. A second set of models, constructed and smeared in the same 
manner as the first, was tested in order to study the sensitivity of the response to different 
keel sizes. The only differences between the two sets are the size of the keel and the size 
of the lumped masses required to maintain the daft of the box structure at one-foot. Keel 
cross section of the second model set is the same as that of the stiffeners, 0.125-in wide 
by 2-in high. Fig. 10 shows a cross section of the small keel model. In both model sets, 


the stiffeners are smeared in a locally orthotropic manner, with an equivalent elastic 
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modulus calculated in the longitudinal and lateral directions. For the second model set, 


the lumped masses were 0.1656-Ibf s*/in’ each. 
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Figure 10. Small Keel Box Model Dimensions. A) Stiffener and Keel Cross 
Sections B) Model Cross Section 


ie First Reduction Method, Both Model Sets 

In the first box model reduction, all stiffeners on the bottom were smeared into 
the shell plating except for the keel and all stiffeners on the sides were smeared into the 
shell plating except for the center longitudinal stiffeners. The bottom, sides, and ends 
and bulkheads were divided into smaller plates so that all stiffeners on those plates would 
be smeared. The reason for this was to obtain a uniform smearing in the model material 
by avoiding the effects from the moment of inertia of the unsmeared keel and stiffeners. 
The equivalent elastic moduli for the two end compartments were considered first. Due 
to the similar dimensions of the sides and bottom of the end compartments, each side and 
bottom of both end compartments was divided into two 12-in by 36-in plates, Fig. 11(A). 


Each plate had two 36-in longitudinally oriented stiffeners and three 12-in laterally 


ie? 


oriented stiffeners smeared into it. This results in each plate having the same longitudinal 
and lateral equivalent elastic modulus. Thus, the sides and bottom of both end 
compartments had the same equivalent elastic moduli. Similarly, the sides and bottom of 
the center compartment were each divided into two 12-in by 48-in plate sections, Fig. 
11(B). Into each plate were smeared two 48-in longitudinally oriented stiffeners and four 
12-in laterally oriented stiffeners. Since each plate was smeared in a like manner, the 


sides and bottom of the center compartment have the same equivalent elastic moduli. 





Figure 11. Smearing of Box Model Plates for the First Reduction Method. 
A) Bottom and Sides of End Compartment B) Bottom and Sides of Center 
Compartment C) Ends/Bulkheads 


For the ends and bulkheads, all stiffeners were smeared into the shell plating 
except for the two center stiffeners. The equivalent elastic moduli were determined by 
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dividing each end and bulkhead into four 12-in square plates, Fig. 11(C). On each plate, 
two 12-in stiffeners were smeared in both the longitudinal and lateral directions. By 
geometry, the four 12-in square plates had the same equivalent elastic moduli in each 
direction. It then followed that each end and bulkhead had the same equivalent elastic 
modulus. The stiffener smearing for the first reduction of the box structure is shown in 


Figee! 22 


Equivalent mass densities for the box structure were calculated in a similar 
manner by dividing the compartments, ends, and bulkheads into smaller plates. The 
equivalent density of each plate was determined by the product of the density and volume 


of each stiffener on the plate divided by the volume of the plate as described above. This 


Bulkheads 





Figure 12. The First Reduction Method. A) Stiffener Smearing on the Sides and 
Bottom B) Stiffener Smearing on the Ends and Bulkheads 
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value was added to the original plate density, yielding the equivalent plate density. Since 
the only difference between the model sets was the unsmeared keel, the equivalent elastic 
moduli and densities for both first reduction models were be the same. Table 1 lists the 
equivalent elastic moduli and mass density values for the first reduction of the box 
models. Fig. 13 shows stiffener location for the first model reduction. The different 
colors of the figure denote sections of the model with the same elastic modulus and 
density. As a result of smearing the stiffeners, the first reduction models have 87 beam 


elements but still have 386 nodes and 378 quad elements. 










a nd 
Bottom and Sides Bottom and Sides 
Table 1. Equivalent Elastic Moduli and Mass Densities for 
the First Reduction Method, Both Model Sets 


Zs Second Reduction Method, Both Model Sets 

For the second box model reduction, all stiffeners on the model except for the 
keel were smeared into the shell plating. Considering the sides of the two end 
compartments, stiffeners on the entire 24-in by 36-in plate were smeared in order to 
calculate the equivalent elastic moduli, Fig. 14(A). This includes five 36-in 
longitudinally oriented stiffeners and three 24-in laterally oriented stiffeners. The bottom 
of the end compartments was smeared in the same manner as the first reduction model by 
dividing it into two 12-in by 36-in plates, Fig. 14(B). The equivalent elastic moduli of 
the center compartment were determined by considering the sides and the bottom 
separately, Fig. 14(C) and (D), respectively. All stiffeners on the 24-in by 48-in side 
plates, five 48-in longitudinally oriented stiffeners and four 24-in laterally oriented 
stiffeners were smeared. The bottom of the center compartment was smeared like the 
bottom of the end compartments by using the 12-in by 48-in plates in the calculations. 
The equivalent elastic moduli for the ends and bulkheads were calculated using the 24-in 


square plate section for smearing with five 24-in stiffeners in the longitudinal and lateral 
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Figure 13. Stiffener Placement for the First Model Reduction 
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Figure 14. Smearing of the Box Model Plates for the Second Reduction Method. 
A) Sides of the End Compartments B) Bottom of the End Compartments 
B) Sides of the Center Compartment D) Bottom of the End Compartment 
E,) Ends/Bulkheads 
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directions, Fig. 14(E). The stiffener smearing for the second reduction of the box 


structure 1s shown in Fig. 15 and stiffener location is shown in Fig. 16. 


The equivalent mass densities were determined in a similar manner by dividing 
the compartments, ends, and bulkheads into smaller plates. As with the first reduction, 
the equivalent density of each plate was determined by adding the original plate density 
to the product of the density and volume of each stiffener on the plate divided by the 
volume of the plate. Like the first reduction, equivalent elastic moduli and densities for 
both second reduction models will be the same. The different colors of the figure denote 
sections of the model with the same elastic modulus and density. Table 2 lists the 
equivalent elastic moduli and mass density values for the second reduction of the box 
models. As a result of the second reduction method, these models had 13 beam elements 


but still had 386 nodes and 378 quad elements. 
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Figure 15. The Second Reduction Method. A) Smearing on the Bottom 
B) Smearing on the Sides C) Smearing on the Ends/Bulkheads 
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Figure 16. Stiffener Location for the Second Model Reduction 
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Compartment Bottom Compartment Sides 
Table 2. Equivalent Elastic Moduli and Mass Densities for 
the Second Reduction Method, Both Model Sets 
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IV. TEST DESCRIPTION AND RESULTS 


A. TEST DESCRIPTION 

All six box structures were subjected to the same shock test. A 20-lb spherical 
charge of HBX-1 was detonated amidships on the starboard (right) side abeam of the 
model at a horizontal range of 8.37 ft and depth of 15.5 ft for a resulting range of 16.75 ft 


as seen in Fig. 17. 
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Figure 17. Shock Test Geometry. A) Top View 
B) Cross Section View 
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B. RESULTS 

Several different nodes throughout the after half of the box model were chosen in 
order to study the overall response trends of the keel, sides, and bulkheads to the shock 
test, as seen in Fig. 18. Seven nodes were selected for examination and represent the 
same location on each of the six models tested. Since the explosive charge was detonated 
at the longitudinal center of the model, the nodal responses of the forward half of the 


structure will be the same as those of the after half due to symmetry of the model. 


A time step of | x 10° seconds was used in the USA/LS-DYNA analysis. This 
time step, based on a characteristic length taken from the longest side of the quadrilateral 
elements, was computed using LS-DYNA. The node responses were recorded at an 
interval of 1 x 10° seconds for 3,000 time steps, representing a total time of 30 
milliseconds. The velocity profiles are shown out to 15 milliseconds in order to study the 


initial velocity responses of the models clearly. 


1. Large Keel Model Results 

Velocity responses of the large keel models are displayed in Figs. 19 — 25. The 
velocity profile for node 55 (see Fig. 18 for location) 1s shown in Fig. 19. Peak 
athwartships velocity of the full-unsmeared model, 255.42 in/sec, was much greater than 
the first and second reduced models, 68.46 in/sec and 86.94 in/sec, respectively. The 
peak resultant velocity responses of the reduced models, 194.47 in/sec for the first 
reduction and 206.11 for the second reduction, were significantly less than that of the full 
model, 264.97 in/sec. It would be expected that the reduced model velocity profiles 
would have greater peak values than those of the full model since the reduced models 
have fewer stiffeners than the full model and are thus not as rigid. In this particular case, 
though, the opposite result has occurred. As can also be seen, velocities of the reduced 


models did not oscillate as much as compared to the full model. 


The athwartships velocity profile for node 66 (see Fig. 18 for location), as seen in 
Fig. 20, shows peaks that were closer in value between the full model and reduced 


models than was the case for node 55, but were still considerably different. The full 
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Figure 18. Node Locations for the Shock Test. A) Keel Nodes 


B) Side and Bulkhead Nodes 








model had a peak velocity of 393.09 in/sec while the first and second reductions had 
453.07 in/sec and 444.58 in/sec, respectively. Peak values differences for the resultant 
velocity were slightly greater with the full model at a peak velocity of 393.38 in/sec and 
the first and second reductions at 480.59 in/sec and 472.39 in/sec, respectively. Both the 
full and reduced models responded in the expected manner. The reduced models had 
greater peak velocities than the full model as a result of having a fewer number of 
stiffeners through out the model, making it less stiff than the full model. Like node 55, 
these responses of the reduced models did not oscillate as much as those of the full 
model. It can be seen that the athwartships and resultant velocities began to converge 
near the same values for all three models, 9 in/sec for the athwartships velocity and 38 
in/sec for the resultant velocity. Also of note, node 66 of the full model had greater peak 
velocities than node 55, but oscillated less freely than node 55. This was due to the fact 
that node 66 has more structural support from the nearby bottom and after bulkhead than 


node 55, which 1s located closer to the top of the structure’s side, as seen in Fig. 18. 


The velocity profiles of keel node 86, shown in Fig. 21, had significantly different 
peak values between the full model and reduced models. Peak vertical velocity of the full 
model was 88.14 in/sec while peak velocity of the first and second reductions were 
209.88 in/sec and 213.10 in/sec, respectively. It would be expected that the keel] node of 
the full model would exhibit a positive initial vertical velocity. However, the vertical 
velocity for node 86 became negative initially, then became positive. Thus, in this 
particular testing case, the keel of the full model deformed in the negative direction 
before deforming in the positive direction. Both reduced models displayed a similar 
behavior that was much less pronounced. The large difference in peak resultant 
velocities is also seen in the velocity profile of node 86. The full model had a peak 
resultant velocity of 110.80 in/sec while the first and second reduction models had peak 
velocities of 210.17 in/sec and 213.32 in/sec, respectively. Again, the models responded 
in the expected manner with the reduced models having greater peak velocities than the 
full model. It can be noted that the full and reduction model velocities quickly converged 
to similar values of approximately 32 in/sec for the vertical velocity and 36 in/sec for the 


resultant velocity. 
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The after bulkhead velocity profiles for node 140 are displayed in Fig. 22. The 
peak velocity for the full model in the longitudinal direction was much greater than the 
peak velocities of the reduced models. In addition, the full model velocity response 
oscillated with greater amplitude than the reduced models. Peak longitudinal velocity for 
the full model was 88.55 in/sec while peak velocities for the first and second reduced 
models were 5.34 in/sec and 9.20 in/sec, respectively. In contrast, the peak resultant 
velocities for the reduced models, 267.56 in/sec for the first reduction and 266.26 in/sec 
for the second reduction, were much greater than that of the full model, 122.46 in/sec. 
The resultant velocity profiles show the expected results with the reduced model peak 
responses greater than that of the full model. However, the longitudinal responses did not 
show the expected behavior. The resultant velocities of both reduced models quickly 


converged near similar values of 35 in/sec. 


Node 157, located on the starboard side of the center compartment as seen in Fig. 
18, experienced higher peak velocities in the reduced models than in the full model, as 
Shown in Fig. 23. The first reduction model had a peak athwartships velocity of 608.40 
in/sec and the second reduction model had a peak athwartships velocity of 593.49 in/sec. 
Peak athwartships velocity for the full model was 344.99 in/sec. Resultant velocity 
responses were similar. The first reduction model had a peak resultant velocity of 626.67 
in/sec and the second reduction had a peak resultant velocity of 612.35 in/sec while the 
full model had a peak resultant velocity of 346.19 in/sec. Here, the models responded in 
the expected manner with the reduced models having greater peak velocities than the full 
model. The full model and both reduced models converged to approximately the same 
athwartships and resultant velocities, 17 in/sec and 35 in/sec, respectively. As can be 


seen, velocities in the reduced models did not oscillate as did those in the full model. 


Node 166 on the starboard side, center compartment showed behavior similar to 
that of node 55. The velocity responses are displayed in Fig. 24. The full model peak 
athwartships and resultant velocities were greater than both reduced models. In addition, 
the full model displayed more oscillation than did the reduced models. The peak 
athwartship velocity for the full model was 333.47 in/sec while the peak athwartships 


velocities for the first and second reduction models were 222.35 in/sec and 249.86 in/sec, 
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respectively. The resultant velocity for the full model was 333.72 in/sec and the resultant 
velocities for the first and second reduction models were 248.11 in/sec and 272.39 in/sec, 
respectively. For node 166, the reduced models did not respond in the expected manner. 
As is seen, the reduced model peak velocities were less than the peak value for the full 
model. Again, velocities for all models in both profiles converged near similar values, 12 
in/sec for the athwartships velocity and 35 in/sec for the resultant velocity. As with the 
comparison between node 55 and 66, node 157 of the full model had greater peak 
velocities than node 166, but oscillated less freely than node 166. This was due to the 
fact that node 157 has more structural support from the nearby bottom and after bulkhead 


than node 166, as was the case with nodes 55 and 66. 


The velocity profiles for keel node 199, shown in Fig. 25, show large differences 
in the peak velocities between the full model and reduced models. The peak vertical 
velocities for the first and second reduction models were, respectively, 321.22 in/sec and 
321.94 in/sec. The peak vertical velocity for the full model was 126.0] in/sec. Like node 
86, the vertical velocity of the full model initially became negative before becoming 
positive. The reason is due to the initia} downward motion of the keel as described for 
node 86. A large difference between the responses is seen in the resultant velocity profile 
as well. The first and second reduction models had peak resultant velocities of 326.65 
in/sec and 327.53 in/sec, respectively, while the full model had a peak resultant velocity 
of 151.19 in/sec. It can be seen that the reduced model peak velocities were greater than 
that of the full model, as expected. Velocities for both profiles converged near similar 
values of 30 in/sec for the vertical velocity and 35 in/sec for the resultant velocity. 


Tables 3 and 4 list the peak velocity responses for the large keel model set. 


Das Small Keel Model Results 
The small] keel mode] set was tested to determine the sensitivity of the velocity 
responses to different keel sizes. Overall, velocity response profiles for the small keel 


models, shown in Figs. 26 — 32, were similar to those of the large kee] models. 
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The peak athwartship velocity for node 55 of the full model, 254.16 in/sec as seen 
in Fig. 26, was much greater than those of the first and second reduction models, 69.99 
in/sec, and 89.08 in/sec, respectively. The resultant peak velocities were somewhat 
closer in value with the full model at 260.08 in/sec and the first and second reduction 
models at 226.15 in/sec and 234.43 in/sec, respectively. As with the large keel model set, 
the reduced models have smaller peak velocities and show much less oscillation than the 


full model. 


The velocity profiles of node 66 are shown in Fig. 27. Results of the small keel 
model set were again similar to those of the large keel model. In the athwartships 
direction, the first reduction had a peak velocity of 454.58 in/sec and the second 
reduction had a peak velocity of 446.10 in/sec while the full model had a peak velocity of 
392.45 in/sec. The resultant velocity profile shows a slight increase in the difference 
between peak velocity values. Peak velocities for the first and second reduction models 
were, respectively, 482.86 in/sec and 474.57 in/sec, and the full model had a peak 
resultant velocity of 392.69 in/sec. Peak velocity behavior is the same as that for node 66 
on the large keel model. The full and reduced model velocities began to converge to 
similar velocity values, approximately 10 in/sec for the vertical component and 33 in/sec 
for the resultant velocity. Comparing the response of node 55 to node 66, node 66 
oscillated much less than node 55 as a result of more structural support from the bottom 


and after bulkhead nearby. This same observation was made for the large keel model set. 


Fig. 28 displays the velocity profiles for keel node 86. Overall, the peak 
responses were less than those of the large keel model set. A significant difference in 
peak vertical velocities still exists between the full and reduced models, though. The full 
model had a peak value of 76.88 in/sec while the first and second reduction models had 
values of 131.19 in/sec and 131.07 in/sec, respectively. Another similarity between the 
large and small keel model sets was the initial negative velocity values for node 86 on all 
three models. Peak resultant velocities are in somewhat closer agreement with the full 
model having a peak of 102.50 in/sec, and the first and second reduction models having 
peaks of 132.31 in/sec and 132.37 in/sec, respectively. As expected, the reduced models 
had greater peak velocities than the full model. Full and reduced model velocities in both 
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vertical and resultant velocity profiles converged to similar values, approximately 32 
in/sec for the vertical velocities and approximately 35 in/sec for the resultant velocities. 


These values are very similar to the values of the large keel set. 


The velocity profiles of the after bulkhead node, node 140, shown in Fig. 29, 
show large oscillations and similar peak velocity responses as was the case in the large 
keel models. However, peak velocities in the small kee] models were more varied than in 
the large kee] model set. Peak velocity for the full model in the longitudinal direction 
was 48.80 in/sec, which ts less than that of the large keel model. The fore and aft peak 
velocity for the first reduction model was 26.19 in/sec, a value greater than the large keel 
model (5.342 1n/sec). The second reduction model displayed a peak resultant velocity of 
9.39 in/sec, a value almost identical to that of the large keel model (9.21 in/sec). The 
reduced models show a greater peak resultant velocity than the full model. The peak 
resultant velocities for the first and second reduction models were 225.98 in/sec and 
226.54 in/sec, respectively while the full model peaked at a velocity of 131.33 in/sec. 
The resultant velocities converged near similar values of 35 in/sec. As compared to the 
large keel model set, the full model peak resultant velocity for the small keel was more 
than that for the large keel, while both reduced small keel model resultant velocities were 


somewhat smaller than those of the large keel models. 


The athwartships and resultant velocity profiles for starboard side node 157, as 
shown in Fig. 30, are similar to those of the large keel model set. Peak athwartships 
velocities for the reduced models were 604.76 in/sec for the fist reduction model and 
589.75 in/sec for the second reduction model. The peak athwartships velocity for the full 
model was 350.13 in/sec. Resultant velocity responses show a peak for the first reduction 
model of 623.05 in/sec and a peak for the second reduction model of 608.67 i1n/sec. The 
full model reached a peak velocity of 351.27 1n/sec. Here, the models responded in the 
expected manner, with the reduced models having greater peak velocities than the full 
model. Velocity values converged near 14 in/sec for the athwartships velocity and near 
35 in/sec for the resultant velocity. Similar to the large keel model set, the velocity 


responses of the reduced models did not oscillate as much as those of the full model. 
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Velocity profiles for node 166, shown in Fig. 31, are also very similar to the large 
keel responses. The full model peak athwartships and resultant velocities were greater 
than those of both reduced models. Additionally, the full model showed more oscillation 
than did the reduced models. The peak athwartships velocity for the full model was 
331.91 in/sec while the peak athwartships velocities for the first and second reduction 
models were 222.37 in/sec and 249.97 in/sec, respectively. Resultant velocity profiles 
show the full model peak velocity of 332.10 in/sec, while those of the first and second 
reduction models were 248.55 in/sec and 272.84 in/sec, respectively. Similar to the large 
keel model set, the reduced models did not respond as expected, having peak velocities 
less than that of the full model. The velocities for all models again converged near 
similar values, 12 in/sec for the athwartships velocity and 37 in/sec for the resultant 
velocity. These values are very close to those found on the large Keel models. As 
explained for the large keel model, node 166 of the full model oscillated more freely than 


did node 157. 


Peak velocity values for node 199 on the small keel model, shown in Fig. 32, 
were greater than those of the large keel model set. The peak vertical components of the 
first and second reduction models were 478.38 in/sec and 479.01 in/sec, respectively, 
while the peak vertical velocity of the full model was 176.16 in/sec. In comparing this 
full model node to the same in the full model with the large keel, the initial vertical 
velocity again became negative before becoming positive. The resultant velocities for the 
small keel models were greater than those of the large keel models. The peak resultant 
velocities for the first and second reduction models were 482.94 in/sec and 483.65 in/sec, 
respectively, while the peak resultant velocity for the full model was 212.23 in/sec. As 
with the large keel model, the reduced model peak velocities were much greater than the 
peak velocity of the full model. Velocities for both profiles converged near similar 
values of 33 in/sec for the vertical velocity and 35 in/sec for the resultant velocity. These 
results are, again, similar to the large keel model set. Tables 5 and 6 list the peak velocity 


responses for the small keel model set. 


The results obtained from smearing the large keel model were varied. For the 


nodes chosen on this model set, peak velocities between the full model and the reduction 
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models were significantly different. Nodes 55, 140, and 166 did not show the expected 
results with regards to the reduced models. These three nodes on the reduced models had 
peak velocities that were less than the peak velocity of the full model. In addition, the 
side nodes of the full model, nodes 55, 66, 157, and 166 as seen in Fig. 18, oscillated with 
greater amplitude than those nodes on the reduced models. This behavior was also 
observed in the response of node 140. Keel node responses of the reduced models (nodes 
86 and 199) show major differences in velocity profiles as well. As can be seen from 
these velocity profiles, smearing seems to have a dampening effect on the velocity 
responses in this particular case. One interesting note is that in most of the responses 
discussed, the first and second reduction models had very similar peak velocity values 


and profile shapes. 


As for keel sensitivity, 1t was noted that most response curves for the small keel 
models were very similar in peak value and shape as those of the large keel models. Like 
the large keel model set, nodes 55, 140, and 166 did not show the expected results of the 
reduced models, having greater peak velocities than that of the full model. There were 
significant differences between the keel and bulkhead responses of the large keel and 
small kee] models. The velocity response of node 86 on the small keel model was less 
sensitive than the large keel model, having smaller peak velocity values than those on the 
large keel model. Keel node 199 showed different behavior. Node 199 on the small keel 
model was more sensitive than the large keel model and displayed greater peak values 
than those on the large keel model. Node 140 displayed different sensitivity behaviors. 
The full model showed a less sensitive longitudinal velocity response while the first 
reduction showed greater sensitivity. The sensitivity of node 140 on the second reduction 
model was unchanged. Conversely, the full model showed more sensitivity in the 
resultant velocity while both reduction models were less sensitive. Like the large keel 
set, most responses of both reduction models were similar in peak velocity value and 


shape. 
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Figure 19. Large Keel Model, Node 55 Velocity Profiles 
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Figure 20. Large Keel Model, Node 66 Velocity Profiles 
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Figure 23. Large Keel Model, Node 157 Velocity Profiles 
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Figure 24. Large Keel Model, Node 166 Velocity Profiles 
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Figure 25. Large Keel Model, Node 199 Velocity Profiles 
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Figure 26. Small Keel Model, Node 55 Velocity Profiles 


45 








Velocity Cin/sec) 


Velocity Cin/sec) 


Node 66 Velocity 


Athwartships 
SOO 


 teestton Model 
 teestton st Reduction 
a 20 eS Reduction 





O 0.003 0.006 0.009 0.012 0.015 


Time (sec) 


Node 66 Resultant Velocity 


—- Full Model 
lst Reduction 


300 


200 


100 





Time (sec) 


Figure 27. Small Keel Model, Node 66 Velocity Profiles 
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Figure 28. Small Keel Model, Node 86 Velocity Profiles 
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Figure 29. Small Keel Model, Node 140 Velocity Profiles 
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Figure 30. Small Keel Model, Node 157 Velocity Profiles 
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Figure 31. Small Keel Model, Node 166 Velocity Profiles 
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Figure 32. Small Keel Model, Node 199 Velocity Profiles 
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Node 140 Fore and Aft (in/sec) 48.80 26.19 
Node 157 Athwartships (in/sec) 350013 604.76 589.75 
Node 166 Athwartships (in/sec) Sey POON 249.97 
Node 199 Vertical (in/sec) 176.16 478.38 479.01 
Table 5. Peak Directional Velocity Values for the Small Keel Model Set 
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V. CONCLUSIONS AND RECOMMENDATIONS 


This thesis examined the effects of smearing beam stiffeners into plates on a ship- 
like box structure that was subjected to an underwater explosion. As the velocity profiles 
showed, responses for the large keel model were significantly different between the full 
model and the reduced models in both peak values and oscillation amplitude. Reducing 
the keel size produced results that were similar to the larger keel case, with the exception 
of the keel and after bulkhead nodes. These nodes did show sensitivities in velocity 


response but did not follow any particular trend. 


It can be concluded, then, that smearing was not effective in this particular case. 
The finite element mesh of this box structure was too coarse to produce results that were 
more comparable between the full model and reduced models. It is considered likely that 
a more finely meshed box model will display more reliable response results. Computer 
simulation of surface ship-shock trials is becoming more attractive in view of the rising 
costs of actual live fire tests and improved smearing techniques will undoubtedly provide 


more acceptable response results for the underwater shock problem. 


Recommended continuing studies in the following areas: 
1. Increase the mesh complexity of the finite element box model. 
2. Increase the geometric complexity of the box structure to study the smearing 
effects of more complex models. 
3. Model the surrounding fluid in order to include the effects of cavitation on the 


smeared models. 
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APPENDIX. USA INPUT DECKS FOR BOX MODEL TEST 


The following programs are the FLUMAS, AUGMAT, and TIMINT input decks 


for the USA code, [Ref. 11]. 
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